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ABSTRACT 

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H, 16 0, H 2 18 0, and H 2 17 - most for the first time 
- in both emission and absorption toward Orion KL using HerschelfHWl. We show that a simple fit, constrained to match the known emission 
and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured 
H 2 18 line fluxes, which are less affected by line opacity than their H 2 16 counterparts, and an escape probability method, the column densities of 
H 2 ls O associated with each emission component are derived. We infer total water abundances of 7.4xl0~ 5 , 1.0xl0~ 5 , and 1.6xlO~ 5 for the plateau, 
hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water 
abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than 
water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from 
(and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources. 

Key words. ISM: abundances — ISM: molecules 



1. Introduction 

During its 6-year mission, the Submillimeter Wave Astronomy 
Satellite (SWAS) surveyed more than 300 galactic sources and 
more than 6800 lines-of-sight (Melnick et al. 2000a), yet none 
produced stronger water emission than the line of sight toward 
Orion-KL. The source of this emission was attributed primarily 
to the chemistry and excitation accompanying the exceptionally 
powerful outflows emanating from the BN/KL region (Harwit 
et al. 1998; Wright et al. 2000; Melnick et al. 2000b; Cernicharo 
etal. 2006; Lerate etal. 2006); however, many sources possessing 
physical conditions favorable to the production of strong water 
emission - e.g. high densities and temperatures - are known to 
exist close to KL and could very likely be significant contribu- 
tors to the water emission detected by ISO, SWAS, and Odin. 
Unfortunately, with access to only the ground-state lio-loi 
transition of ortho-H 2 6 Opland H 2 ls O, even the velocity -resolved 
SWAS and Odin measurements were limited in what could be 
inferred about the various components giving rise to the strong 
water emission. 

The availability of the HerschelfHWl instrument (deGraauw 
et al. 2010) with its extended frequency coverage and higher an- 
gular resolution, now permits a more detailed examination of the 
conditions responsible for the water emission toward Orion-KL. 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 

1 Also referred to simply as H2O 



Here we report the detection of 21 H 2 0, 15 H 2 ls O, and 12 rL/'O 
velocity-resolved lines toward this source obtained as part of the 
HEXOS program (Bergin et al. 2010). 

In this paper, we present an analysis of the sources of the 
water emission based upon the lower-opacity lines of H 2 18 0. We 
also show that the approach taken in this analysis holds great 
promise when applied to the H2O and H 2 n O lines, which will be 
pursued in a future paper. 



2. Observations and results 

The HIFI observations presented here were carried out in March 
and April 2010 using the spectral scan dual beam switch (DBS) 
mode pointed towards Orion-KL 072000 - 5 /! 35'"14.3 i and 
<5/20oo = -5°22'33.7". All observations were obtained with a 
beamsize of ~ (22 / vthz)" and reference beams approximately 
3' east and west, which is roughly orthogonal to the orienta- 
tion of the Orion molecular ridge (e.g., Ungerechts et al. 1997). 
However, water emission is extended in Orion (Snell et al. 2000) 
and the reference beam may contain some contamination from 
a narrow (Av~ 3-5 km s _1 ) component centered at ~9 km 
s . We utilized the wide band spectrometer providing a spectral 
resolution of 1. 1 MHz over a 4 GHz IF bandwidth. The data pre- 
sented here are from a range of HIFI bands obtained as part of 
the HEXOS program. These data were reduced and converted to 
single side band as described by (Bergin et al. 2010), with ad- 
ditional analysis performed at the CfA. In our study, we adopt a 
uniform main beam efficiency of 70%. 
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Because of flux differences between the H- and the V- 
polarizations, which are most likely due to the known pointing 
offset between the two beams, we use only the H-polarization 
data for our analysis. The spectra for all H2O, H 2 17 0, and H 2 18 
lines were extracted from the more extended HEXOS spec- 
tral scan data using the JPL Spectral Line Catalog (Pickett 
et al. 1998) for identification. Finally, the continuum offset ap- 
propriate to each line was determined directly from emission- 
free spectral regions near each line. 

Figures 1 and 2 show the spectra of H 2 8 and H2O plus 
H 2 17 0, respectively. These spectra span a broad range of exci- 
tation conditions, ranging in upper-level energies between 53 K 
and more than 1000 K. All spectra have been examined for se- 
vere blending using the CLASS-Weeds tool (Maret et al. 2010), 
the JPL Spectral Line Catalog, or visual evidence of non-smooth 
water line wings. Blended lines were excluded from the follow- 
ing analysis. 



3. Analysis 

The goals of the present effort are twofold: (1) isolate the compo- 
nents giving rise to the water emission we detect; and, (2) model 
these components in a way that best reproduces the measured 
line fluxes and profiles. To do this, we focus here on the observed 
H 2 ls O lines. These lines have been detected over a broad range 
of excitation conditions with high signal-to-noise ratios and are 
much less affected by optical depth effects than their H 2 6 coun- 
terparts, making the analysis more straightforward. In addition, 
the 16 0: 18 ratio is well known (i.e., ~ 500) and not believed to 
vary significantly between sources, making the conversion from 
inferred H 2 s O abundance to H 2 16 abundance robust. 

Step 1 - isolate the components: The lines exhibit complex 
profiles which we attribute to a combination of emission and ab- 
sorption components along the line of sight. To isolate what we 
believe are the three predominant emission components within 
the HIFI beams - namely the plateau molecular outflow, the hot 
core, and an extended region of gas composed of the compact 
ridge plus the warmer, denser portion of the extended ridge near 
KL (cf. Blake et al. 1987) - we adopt a line-fitting strategy that 
fixes the well-established characteristics of these regions, such as 
their v L sr, and, in some cases, the typical line width, and leaves 
as free fitting parameters such quantities as the line strengths. 

In addition to the three emission components, we include 
the effects of absorption by foreground material in two dis- 
tinct kinematic components: a narrow component near 7 km s -1 , 
and a broad component centered at an LSR velocity of -5.1 
km s -1 . While the presence of these absorption components is 
clearly required to fit the observed water line profiles, particu- 
larly in the case of low-lying transitions of H 2 16 0, the existence 
of foreground absorbing material at these velocities has been 
independently confirmed by HIFI observations of HF (Phillips 
et al. 2010), OH + and H 2 + (Gupta et al. 2010), as well as 
CRIRES observations of the fundamental CO vibrational band 
(Beuther et al. 2010). The narrow component arises in quies- 
cent gas, while the broad, blueshifted component represents out- 
flowing material, presumably associated primarily with the Low 
Velocity Flow (Genzel & Stutzki 1989). For the lower-lying tran- 
sitions, these absorption components account for pronounced 
asymmetries in the line shapes, as well as the absorption fea- 
ture close to the systemic source velocity (although we note here 
that narrow line emission in the reference beam is potentially a 
contributor to this absorption feature observed in the very low- 
est transitions.) Even in the case of H 2 18 0, transitions to the 



Table 1. Fixed and varied parameters in water-line fits 



Source 



Peak T* A 
(K) 



V LSR 

(km s" 1 ) 



FWHM 

(km s _1 ) 



Plateau Varied +6.9 Varied 

Hot Core Varied +5.2 10.0 

Extended warm gas Varied +8.25 2 - 8 1 

Narrow absorption Varied +6.88 6.70 

Broad absorption Varied -5.1 30.0 

' The FWHM was constrained to vary between only 2 and 8 km s _1 . 

ground states of ortho- or para-H 2 8 (i.e., 2^ - lot, Lo _ loi 
and 1 1 1 - Ooo) are affected by foreground absorption. Indeed, in 
the In -Ooo and 2i2 - lot H 2 18 transitions, where the contin- 
uum brightness temperature is greatest, the blueshifted absorp- 
tion feature can cause the observed antenna temperature to dip 
below the continuum level. 

Thus, fits to all lines were made using the expression: 

Fitted Line (T A *) = (Continuum Offset + G p i at + Gh c + G ewg ) 

x Exp[-(G na + G ba )], (1) 

where G p i at , Gh c , G eW g, G na and Gb a are Gaussian components 
representing the plateau, hot core, extended warm gas, narrow 
absorbing feature, and broad absorbing feature, respectively. 
Table 1 provides the fit parameters fixed by previous measure- 
ments and those that were allowed to vary, unconstrained, in or- 
der to obtain the best fit to the line profiles. 

Step 2 - model the H 2 s O emission components: The results 
of Step 1 are a set of best-fit integrated intensities for each com- 
ponent and transition, including the absorption features, that sum 
to reproduce the line flux and profile for each ortho- and para 
H 2 18 line. In this paper, we focus on the emission components 
only; analysis of the physical conditions associated with the ab- 
sorption components will be undertaken following the results of 
a soon-to-be-completed water map toward Orion-KL. To assess 
how the H 2 18 line strengths constrain the water abundance in 
each component, the equilibrium level populations of all H2O 
ortho and para rotational levels of the ground vibrational state 
with energies E/k up to 2000 K have been calculated using an 
escape probability method that includes the necessary effects 
of radiative excitations due to dust emission embedded within 
each component. It is assumed that the water molecules see An 
steradians of dust emission from within each component. The 
velocity gradient for each transition is assumed to be equal to 
Av n(H2)/V(H2), where the line width, Av, for each line for each 
component is taken from the best fit in Step 1, and n(H2) and 
V(H2) are the volume and column densities of H2, respectively. 
The rate coefficients for collisions between ortho- and para-H2 
and ortho- and para-H20 calculated by Faure et al. (2007) are 
used, and the H2 ortho-to-para ratio is assumed to be the LTE 
value at the gas temperature of each component. Finally, the cal- 
culations incorporate the beam size and aperture efficiency ap- 
propriate to each transition. 

More than 90% of the presently observed H2O total line flux 
(and > 98% of the H 2 18 and H 2 n O total line flux) lies in transi- 
tions with £upper < 600 K. Thus, we focus our modeling efforts 
primarily on reproducing the flux and profiles for these transi- 
tions. The H2 density, gas and dust temperatures, source size, 
and ortho- and para-H 2 18 column densities were varied to best 
match the inferred line fluxes for each emission component. The 
values yielding the best fit to the data are provided in Table 2. 
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Table 2. Best- fit radiative transfer model parameters for Orion H 2 ls O emission components 

Av r gas »(H 2 ) r dust t A?(ortho/para-H 2 ls O) Total Inferred 

Source (kms -1 ) (K) (cm -3 ) (K) iOuvct , (cm -2 ) H 2 Abundance * 

Extended warm gas 2-8 75 2xl0 6 30 20" 7.4xl0 15 (ortho)/ 2.0xl0 15 (para) 1.6xl0- 5 

Plateau 20-34 188 2xl0 6 113.6 24" 1.2xl0 16 (ortho)/ 2.8xl0 15 (para) 7.4xl0~ 5 

Hotcore 10 150 lxlO 7 180 5" 1.5xl0 16 (ortho) /5.0xl0 15 (para) LOxlO -5 

t Greybody fit to the Orion continuum of the form: B f (T dust ) x (0.0233 y) - 486 , where v is wavelength in wavenumbers. + Assumes 16 0/ 18 = 500 
and NQI2) = 3xl0 23 , lxlO 23 , and lxlO 24 cm -2 in a 30" beam for the extended ridge, plateau, and hot core, respectively (Blake et al. 1987). 



The line profiles resulting from the radiative transfer model cal- 
culations for the emission components and Step 1 line-fits to the 
absorption components are shown as the red curves superposed 
on the observed spectra in Fig. 1. The models summarized in 
Table 2 provide a remarkably good match to the data, though 
the deviation between the models and the observed spectra for 
the higher-energy H 2 I8 transitions clearly illustrates the short- 
comings of single-value models for each component as small 
amounts of hotter gas are not accounted for. 

The physical conditions summarized in Table 2 have also 
been used to model the H 2 16 lines with £uppei <400 K. To do so, 
the column densities of ortho- and para-H 2 16 are assumed to be 
500 times greater than those of H 2 ls O, the line fluxes calculated, 
and then applied using the best-fit H 2 16 plateau line widths de- 
termined using Eqn. 1 . For the hot core and extended warm gas 
region, the H 2 16 widths were assumed to be twice those of the 
H 2 s O, and the absorption components are unchanged. The re- 
sults of this simple approach are shown as the superposed red 
curves on the relevant H 2 16 spectra in Fig. 2. The potential for 
a more careful analysis of the H 2 16 and H 2 n O lines is illus- 
trated by how well the constrained fits match the other line pro- 
files, shown as the superposed brown curves in Fig. 2. A more 
detailed model will be presented in a future paper. 



4. Discussion 

Modeling of the rich spectrum of H 2 ls O lines toward Orion- 
KL reveals several things. First, the relatively high H2O abun- 
dance associated with the plateau is consistent with elevated wa- 
ter abundances measured previously toward KL (cf. Cernicharo 
et al. 2006) as well as toward a number of other molecular out- 
flows (cf. Franklin et al. 2008). This is most likely the result of 
a combination of H20-ice sublimated and sputtered from grain 
surfaces and H2O formed efficiently in the gas phase via neutral- 
neutral reactions favored in hotter portions of the plateau. The 
inferred water abundance for the plateau given in Table 2 is 
less than that cited in some larger-beamsize studies (e.g., Harwit 
et al. 1998, Melnick et al. 2000b), and may be due to the ex- 
clusion of more extended regions where the outflows encounter 
the surrounding quiescent material (cf. Genzel & Stutzki 1989). 
These shock-heated regions, which are particularly prominent in 
H2 emission, can subject virtually all of the affected gas to tem- 
peratures in excess 1000 K, thus facilitating the neutral-neutral 
reactions that efficiently produce H2O. 

Second, the water abundances inferred for the hot core and 
extended warm gas are more than an order of magnitude greater 
than that inferred toward other quiescent regions (cf. Melnick & 
Bergin 2005). This is likely the result of enhanced sublimation 
of water-ice from, and reduced freeze-out onto, the warm dust 
grains present within both regions. It should be noted that the 
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Fig. 1. H 2 18 lines toward Orion-KL in order of increasing up- 
per level energy. The superposed red curves show the line pro- 
files resulting from our radiative transfer modeling of the emis- 
sion components and line fits to the absorption components. 
The labels in the upper left corner of each plot list the species, 
the transition, the transition rest frequency, and the upper-level 
energy. The vertical dashed line denotes the 9 km s -1 sys- 
temic velocity of the cloud. The 22o-2n, 423-330, 634-541, 
and 633 - 542 spectra are omitted due to blending with other 
lines or a low signal-to-noise ratio. We also note possible blend- 
ing of the 4 22 - 4i3 line with CH 3 OH (12 5 1 - ll 2 i) and H 2 13 CO 
(18 2 17 - 18ois), both of which lie within 27 km s -1 of the H 2 ls O 
line. 



Melnick et al.: Water in Orion 




H 2 1,„-1 D1 
556.94 GHz 
61 .OK 




H z O 2 0z -1 „ 
987.93 GHz 
100.8 K 




H z O 2„-2„ z 
752.03 GHz 
136.9 K 




H z O 2 20 -2„ 
1228.79 GHz 
195.9 K 




H z O 3,2-2,, 
1153.13 GHz 
249.4 K 




R (km s" 1 ) 




H z O4„-4 04 
1602.22 GHz 
396.4 K 




H z O 3 3 ,-4„ 
1893.69 GHz 
410.4 K 



yifHf 



mW 



H z 4 zz -4„ 
1207.64 GHz 
454.3 K 




H z O 5 3Z -5 ZS 
1867.75 GHz 
732.1 K 





H z O 6 Z4 -6 IB 
1794.79 GHz 
867.3 K 




H z O 6 s4 -7„, 
1880.75 GHz 
933.7 K 





, (km s" 1 ; 





H z "0 1„ -0„ |, 
1107.17 GHz J \ 
53.1 K 






V, 



Jl 



H,"0 2„-1, 
1662.46 GHz 



Hj"0 2„ - 2 
748.46 GHz 
136.6 K 






H 2 "0 4,, 
1604.18 GHz 
395.9 K 






Fig. 2. Left: Same as Fig. 1 , except showing the H 2 16 spectra toward Orion-KL/Hot Core in order of increasing upper level energy. 
The red curves superposed on the H 2 16 spectra with upper-level energies less than 400 K result from our radiative transfer model 
of the emission components and line fits to the absorption components. The brown curves show the best 5-component fit resulting 
from the procedure described in Section 3. The 22i - 2i2, 32i - 3i2, 624 - 7 17, 734 - 725, and 744 - 651 spectra have been omitted due 
to blending with other lines or a low signal-to-noise ratio. Right: H 2 n O spectra. The brown curves superposed on the spectra show 
the best 5-component fit resulting from the procedure described in Section 3. The 2q2 - 1 11 and 221 - 2i2 spectra have been omitted 
due to blending with other lines or a low signal-to-noise ratio. 



gas and dust temperatures inferred for the extended warm gas 
should be viewed as lower limits given the probable presence of 
both water-line and continuum emission in the reference beam. 
Finally, the H, 18 ortho-to-para ratio inferred for all three 
emission components is consistent with a ratio of 3:1. A ratio of 
greater than 3:1 is likely the consequence of the rather simple 
model adopted for each component or residual inaccuracies in 
the water collisional rate coefficients, or both. 
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Appendix 



Table 3. Best-Fit H 2 18 Integrated Line Intensities t 



Transition 



Frequency 
(GHz) 



Upper-Level 

Energy 

(K) 



Plateau 
jTJdv 

(K-km s- 1 ) 



Emission Component 
Hot Core 

JT;dv 

(K-km s- 1 ) 



Extended 
Warm Gas 

JTJdv 
(K-kms" 1 ) 



lu-Ooo H01.70 

Lo-loi 547.68 

2o2-ln 994.68 

2i2- lot 1655.87 

2 U -2 02 745.32 

2 21 - 2n 1633.48 

3i2-3o3 1095.63 

3i2-2 21 1181.39 

3 2 2-3i 3 1894.32 

3 2 i-3i2 1136.70 

4 13 -4 04 1605.96 

4 22 -4 13 tt 1188.86 

5 24 -4 3 i 1003.28 

5 3 2-5 23 1815.85 

6 24 -6 l5 1800.47 



52.9 
60.5 
100.6 
113.6 
136.4 
192.0 
248.7 
248.7 
294.6 
303.3 
395.4 
452.4 
595.9 
727.6 
865.0 



315.4 
199.0 
423.3 
574.9 
276.6 
220.6 
444.9 
178.8 

65.7 
424.4 

46.6 

118.5 

0.0 

22.2 

41.5 



7.6 

8.2 
60.1 

0.0 
66.3 
10.1 
70.2 
50.4 
34.7 
65.3 
38.6 
57.5 

6.9 
32.8 
36.9 



50.8 

34.6 

10.7 

71.5 

17.7 

6.4 

5.2 

1.8 

2.4 

11.8 

0.0 

0.2 

2.0 

10.0 

0.0 



t The line fitting procedure is described in Section 3. 
tt Some flux attributed to this transition may be due to the CH 3 OH (12 5 1 - 1 1 2 1) and H 2 13 CO (I8217 - 18 is) 



transitions, both of which lie within 27 km s 



oftheH 2 18 4 2 2-4 13 line. 



